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Predictor-Feedback Stabilization of 
Multi-Input Nonlinear Systems 

Nikolaos Bekiaris-Liberis and Miroslav Krstic 


Abstract 

We develop a predictor-feedback control design for multi-input nonlinear systems with distinct input delays, 
of arbitrary length, in each individual input channel. Due to the fact that different input signals reach the plant 
at different time instants, the key design challenge, which we resolve, is the construction of the predictors of the 
plant’s state over distinct prediction horizons such that the corresponding input delays are compensated. Global 
asymptotic stability of the closed-loop system is established by utilizing arguments based on Lyapunov functionals 
or estimates on solutions. We specialize our methodology to linear systems for which the predictor-feedback control 
laws are available explicitly and for which global exponential stability is achievable. A detailed example is provided 
dealing with the stabilization of the nonholonomic unicycle, subject to two different input delays affecting the speed 
and turning rate, for the illustration of our methodology. 

I. Introduction 

a) Background and Motivation: Despite the recent outburst in the development of predictor-based 
control laws for nonlinear systems with input delays (511, O, (71, O, (H, (TOl . (TTI . (T3l . (T4l . (TSl . (Tbl . 

uni, m, m, lEi, m, (sb, m, (n, m, m, m, m, m, (n, m, ^e problem of the 

systematic predictor-feedback stabilization of multi-input nonlinear systems with, potentially different, 
in each individual input channel, long input delays, has remained, heretofore, untackled, although the 
problem was solved in the linear case in the early 1980s 0| (see also iHTI l. In this article, we address the 
problem of stabilization of multi-input nonlinear systems with distinct input delays of arbitrary length and 
develop a nonlinear version of the prediction-based control laws developed in (H and recently in (53l . 
(54ll for the compensation of input delays in multi-input linear systems. 

Besides the unavailability of a systematic predictor-feedback design methodology for multi-input non¬ 
linear systems with long input delays, the real motivation for this article comes from applications. 
Such systems serve as models for the dynamics of traffic (221 . (47l . teleoperators Ii24l and robotic 
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manipulators [l2l, [l20ll . motors ll34l . [l50ll . multi-agent systems (H, ifT^ . [|40l, autonomous ground vehieles 
dsn, unmanned aerial vehieles and planar vertieal take-off and landing airerafts EH . fl48l . and the 
human museuloskeletal system in applieations sueh as neuromuseular eleetrieal stimulation E^ . E^ . EH . 
to name only a few. Motivated by the negative effeets of input delays on the stability and performanee 
of sueh eontrol systems, in this artiele we present eontrol designs that aehieve delay eompensation. 

b) Contributions: We introduee a predietor-feedbaek eontrol design for the eompensation of long 
input delays in multi-input nonlinear systems. Sinee eaeh individual input ehannel might induee a different 
delay the predietors of the plant’s state are eonstrueted reeursively starting from the predietor that 
eorresponds to the smallest input delay all the way through to the predietor that eorresponds to the 
largest input delay. Speeifieally, at eaeh step, the predietor, over the predietion horizon that eorresponds 
to the eurrent’s step input delay, is eonstrueted by aetually predieting, over the appropriate predietion 
window, the future values of the predietor eonstrueted at the previous step. 

We eonduet the stability analysis of the elosed-loop system, under the developed predietor-feedbaek 
eontrol law, utilizing two different teehniques-one based on the eonstruetion of a Lyapunov funetional and 
one based on estimates on the solutions of the elosed-loop system. In the former ease, the eonstruetion 
of a Lyapunov funetional is enabled by the introduetion of novel baekstepping transformations of the 
aetuator states, whieh are based on an equivalent, PDE representation of the eonstrueted predietor states. 
In the latter approaeh, we exploit the faet that eaeh delay is eompensated after a finite time. 

We present a detailed example, ineluding numerieal simulations, dealing with the stabilization of a 
nonholonomie robot subjeet to different input delays, in order to highlight the intrieaeies of our design 
and analysis methodologies. We speeialize our results to the ease of linear systems for whieh the predietor- 
feedbaek eontrol laws are obtained explieitly and for whieh global exponential stability is aehievable. 

c) Organization: We start in Seetion with the introduetion to the problem of predietor-feedbaek 
stabilization of multi-input nonlinear systems and develop the predietor-feedbaek eontrol laws. In Seetion 


III we prove global asymptotie stability of the elosed-loop system under predietor-feedbaek by eonstrueting 


a Lyapunov funetional and in Seetion]^ we prove global asymptotie stability using estimates on solutions. 
Seetion |V] is devoted to a detailed example of stabilization of the nonholonomie unieyele subjeet to input 
delays. We speeialize our methodology to the ease of linear systems in Seetion |V^ For the example 
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Fig. 1. Multi-input nonlinear system with distinct input delays. 


Notation: We use the common definition of class /C, /Coo and ICC functions from [[33l. For an n- 
vector, the norm | ■ | denotes the usual Euclidean norm. For a function u : x IR+ —)■ M we denote 

by ||m(/:)||oo its spatial supremum norm, i.e., ||m(/:)||oo = sup 3 ,g[o,_D] For any c > 0, we denote the 

spatially weighted supremum norm of u by II 'w(f)||c,oo = sup 3 ,g[o,D] t)|. For a vector valued function 

p : X ]R+ —)■ M" we use a spatial supremum norm ||p(t)||oo = sup 3 ,g[o^D] \/pi(a;, + .. .pn{x,tY. 

We denote by C^{A;E) the space of functions that take values in E and have continuous derivatives of 
order j on A. 

IF Multi-Input Nonlinear Systems with Distinct Delays 
AND Predictor-Feedback Controe Design 

We consider the following system (see Fig. 

X{t) = f {X{t), Dm)), (1) 

where X G M” is state, f > 0 is time, f/i,..., t/m £ IR are control inputs, Di,..., Dm are (potentially 
distinct) input delays satisfying (without loss of generality) 0 < Di < ■ ■ ■ < Dm, and / : x MD —)■ M" 

is a locally Fipschitz vector field that satisfies /(O, 0,..., 0) = 0. The predictor feedback controllers are 
defined by 

Ui{t) = Ki{Pi{t)), i = (2) 

where Nj : —)■ M, / = are continuously differentiable feedback laws with Ni(0) = 0, 

i = and. Pi are the Dj-time units ahead predictors of X, for all i = Defining 


worked out in detail in Section M we present simulation results in Section VII 
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Dji = Dj — Di, for all i < j <m, the predictors are given by 

Pi(t) = x{t)+f f{p,{e),u,{e),U2{e-D2i),...,uuo-D,^,))de 0) 

J t—Di 

P2{t) = Pi{t)+f f{P2{e),K,{P2{e)),U2{e),u,{e-D,2),...,uuo-Dm2))de (4) 

J t — D21 


Pmit) = Pm-lit)+ fiPU0),K,iPU0)),K2iPmm,---,f^m-liPmm,UU0))deA5) 

*' ^ ^mm — 1 

with initial conditions for the integral equations Q-Q 


Pl(0) = X(0)+/ f{P,{s),U,{s),U2{s-D2l),...,Um{s-Drr,l))ds, -P, < 0 < Q (6) 

J-Di 

P2{0) = Pl(0)+/ f{P2{s),K^{P2{s)),U2{s)Ms-D^2),...,U^{s-D^2))ds, 

J —D21 


-D 21 <0 <0 


(7) 


Pm{9) = Pm_i(0)+/' f {Pm{s),Ki{Pm{s)) ,H2{Pm{s)) . ,K^-i{P„,{s)) ,Um{s))ds, 

^mm — 1 

-Dmm-l <e<0. (8) 

We show that Pi, for all z = 1,... ,m, are the Pj-time units ahead predictors of X by induction. In 
order to better understand the general induction step we provide two initial steps. We show first that Pi 
and P 2 are the Pi- and P 2 -time units ahead predictors of X, respectively. 

Step 1: We perform the change of variables t = 6 + Pi, for all 9 > —Pi, in ([^ and define Pi{9) = 
X{9 + Pi), 9 > —Pi, to arrive at 
dP (9^ 

= f (Pi(0), U,{9), U2{9 - P 21 ), ...,Um{9- Pmi)), for all 0 > -Pi. (9) 

Integrating (|^ from 9 = t — Di to 9 = t and using definition Pi(6*) = X{9 + Pi) we get Q. Integrating 
(|^ from 9 = —Pi to any 9 <0 and using definition Pi (—Pi) = X(0) we get ([^. 

Step 2: Performing the change of variables 6* = s + P 21 , for all s > —D 21 , in Q and defining 
P 2 (s) = Pi(s + P 21 ) = X(s + P 2 ), for all s > -P 21 , we get that 

= / (P 2 (s), «1 (P 2 (s)), U 2 {s), Us{s - P 32 ),..., Um{s - D^ 2 )), for all 5 > -P 21 , (10) 
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where we also used the faet that Ui{s + D21) = ki {Pi{s + D21)), for all s + D21 > 0, and definition 
F 2 (s) = Pi{s + D21). Integrating ( [T^ from s = t — D2i io s = t and using definition P2{s) = Pi{s + D2i), 
for all s > —-D21, we arrive at Q. Integrating ( [T^ from 6 = —D21 to any 6 * < 0 and using definition 
P 2 {-D 2 i) = Pi(0) we get 0. 

Step j: Assume now that the Dj-time units ahead predietor of X, namely Pj, satisfies the following 
ODE in r 

dPjir) 


dr 


= f {Pj{r )),..., Kj_i {Pj{r )), Uj{r), Uj+i{r - Dj+ij ),..., U^{r - D^j )), 


for all r > — (11) 

Performing the ehange of variables r = h + Dj+ij, for all h > —Dj+ij, in ( [TT] ) and defining Pj^i{h) = 
Pj{h + Dj+ij) = X{h + Dj+i), for all h > —Dj+ij, we get that 
dPj+i{h) 


dh 


= f Ki {Pj+i{h)) , Kj {Pj+i{h)) , + , Um{h - Dmj+l)) , 


for all h > —Dj^ij, 


( 12 ) 


where we also used the faet that Uj{h + = kj {Pj{h + Dj^ij)), for all h + > 0, and 

definition Pj+i(h) = Pj{h + Dj+ij). Integrating (12) from h = t — -Dj+ij to h = t and from h = —Dj+ij 
to any h < 0, and using definition Pj+i{h) = Pj{h + Dj^ij), for all h > (whieh implies 

that = Fj(0)), we eonelude that indeed the Dj-time units ahead predietors of X, for all 

i = 1,... ,m, are given by 0-0 with initial eonditions 0-0. 

III. Lyapunov-Based Stability Analysis Under Predictor Feedback 


Assumption 1: The system X = f (X, wi,... , 0 ;^) is strongly forward eomplete with respeet to a; = 

(cUi, . . . , UJ^n) . 

Assumption 2: The system X = f {X,uji + k,i{X), ... ,Um + UmiX)) is input-to-state stable with 
respeet to ce = (cei,..., corn)'^- 

Assumption 3: The systems X = gj (X, Wj+i, .. . , Um), for all j = 1,..., m— 1, with gj (X, uij+i, ..., Um) = 
f (X, (X ),... ,Kj (X), ujj+i ,..., cOrn), ^Te sfiongly forward eomplete with respeet to a; = (wj+i,..., . 

The definitions of strong forward eompleteness and input-to-state stability are those from OTll (see also 
li3l for the definition of standard forward eompleteness whieh differs from strong forward eompleteness 
in that / (0, 0,..., 0) =0) and [[52l . respeetively. 
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Assumption guarantees that for every initial condition and every locally bounded input signal the 
corresponding solution is defined for all f > 0. In particular, the plant does not exhibit finite escape 
before the first feedback control reaches it. This is a natural requirement for achieving global stabi¬ 
lization in the presence of arbitrary large delays affecting the inputs of a system. Assumption can 
be relaxed to only global asymptotic stability of system X = f {X, ki{X),..., Km{X)). Yet, at the 
expense of not having a Lyapunov functional available. Assumption guarantees that after the j-th 
controller “kicks in” and the Dj-\h delay is compensated, and hence, the plant behaves according to 
X = f (X, Ki{X),..., Kj (X), Uj+i{t — Dj^i),... ,Um{t — Dm)), the solutions are also well-defined. In 
particular, the plant does not exhibit finite escape before the j-fl-th feedback control reaches it and after the 
j-th feedback control has already reached the plant. Note that Assumption can be relaxed to strong for¬ 
ward completeness of systems X = gj {X, coj+i ,..., Um) with respect to a; = (ojj+i ,..., ujm)'^, for all j G 
{ri,ri +r 2 ,...,ri + ...+ r„}, where gj {X,u;j+i,... ,Um) = f (X, ki (X ),..., kj (X) ... ,Wm), 

ri denotes the number of delays that are equal to Di, a = 2,... ,u, denotes the number of delays 
that are equal to and u is the number of distinct delays. In particular, when all delays are 

identical. Assumption can be completely removed. 

The stability proof is based on an equivalent representation of plant Q, using transport PDEs for the 
actuator states, and on an equivalent PDE representation of the predictor states Q-Q. We present the 
alternative representations for the plant and the predictor states before stating and proving the main result 
of this section, since the reader might find the alternative formalisms helpful in better digesting the design 
and analysis ideas of our methodology. 

A. Equivalent Representation of the Plant Using Transport PDEs for the Actuator States 


System Q can be written equivalently as 

X{t) = f{X{t),Ui{0,t),...,UmiO,t)) (13) 

dtUi{x,t) = d^Ui{x,t), xe{0,Di), i = l,...,m (14) 

Ui{Di,t) = Ui{t), i = l,...,m. (15) 

To see this note that the solutions to ( [T4] ), ( [T5] ) are given by 

Ui{x,t) = Uft + X - Di), xe[0,Di\, i = l,...,m. (16) 
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X Pi P 2 Pm-1 Pm 


Pi{x) 

P2{x) 


Pm{x) 

Ui{x), . . .,Um{x) 

U2{x), . . . ,Um{x) 


Um{x) 


U\ K\ {Pi) U 2 ^2 {P 2 ) Pm—1 1 {Pm—l) Pm {Pm) 

0 Di D 2 Dm-l Dm 

I-^^^ 


Fig. 2. The Dj-time units ahead predictors of the state X, namely Pj, given in 0-0, and their equivalent representation by the PDF 
states Pj, given in (|17^-(|19^, based on the transport-PDE equivalent of the actuator states defined in ^14[l-l|15^. The control laws Uj are 
defined in in terms of Pj and can be written equivalently as in in terms of pj. 


B. Transport PDE Representation of the Predictor States 

The predictor states Pi{9), for all 9 > —Di, and Pj{9), for all 9 > —Djj^i and j = 2,... ,m, can be 
written equivalently as (see Fig. 

Pi{x,t) = X{t)+ [ f {pi{y,t),Ui{y,t),...,Um{yP))dy, x G [0, T>i] (17) 

nx 

P2{x,t) = pi{Di,t)+ f {P2{y,t),ni{p2{y,t)) ,U2{y,t),... ,Um{y,t))dy, xg[T>i,T> 2 ] (18) 

JDi 


rx 

Pm{x,t) = Pm-l{Dm-l,t) + f {Pm{y, t), {pm{y, t)) , . . . , Km-l {Pm{y, t)) , Um{y, t)) dy, 

d Dm-l 

X G [Dm-l, Dm]- (19) 

We show this by induction. In order to make the presentation of the procedure clearer we present two 
steps before the general step. We first observe that (see Section Pi{9) = X(9 + Di), for all 9 > —Di, 
and Pj{9) = X{9 + Dj), for all 9 > —Djj_i and j = 2, ...,m. The functions Pi, i = 1, ... ,m, satisfy 
the following ODEs in x 

dxPi{x,t) = f {pi{x,t),ui{x,t),... ,Um{x,t)), xe[0,Di\ (20) 

dxP 2 {x,t) = f{p 2 {x,t),ni{p 2 {x,t)),U 2 {x,t),...,Um{x,t)), XG[T>i,T>2] (21) 


dxPm{x,f) = f {pm{x,t),Ki{pm{x,t)) m ■ ■ ,K-m-l{Pm{x,t)) ,Um{x,t)) , Xe[Dm-l,Dm], ( 22 ) 











with initial conditions 


Pi{0,t) = X{t) 

P2{Di,t) = pi{Di,t) 


(23) 

(24) 




(25) 


Step 1: The solution to (20), (23) is 


Pi{x,t) = X{t + x), xe[0,Di]. 


(26) 


In order to show this, first note that (26) satisfies the boundary condition (23). The function X{t + x) 
also satisfies the ODE in x (20) whieh follows from the fact that by ([^ one can conclude that 

X'{t + x) = f {X{t + x), Ui{t + X — Di),..., Um{t + X — Dm)), for all f > 0 and 0 < x < Di. (27) 


The result follows from the uniqueness of solutions to the ODE ([^. Therefore, by defining 

Pi{Di,t) = Pi{t), 


(28) 


and using the faet that pi is a funetion of one variable, namely x + t (whieh follows from (26)), one ean 
eonelude that 


Pi{t + x-Di) = pi{x,t), xe[0,E>i]. 

Performing the change of variables x = 9 + Di — t, for all t — Di < 9 < t, in ([Tt]) and using ([T^, 


(29) 


and (29) we arrive at 


P,{9) = Xit)+ f{P,{s),U,{s),...,Um{s-Dmi))ds, forallt-D,<9<t. (30) 
J t—D\ 


Step 2: Similarly, it ean be shown that 


P 2 {x,t) = X{t + x), xe[E>i,E> 2 ], 


(31) 


is the solution to ( [2T] ), ( |24l ) sinee it satisfies ( |24l ) and sinee it satisfies the ODE in x ( |2T] ) whieh follows 
from the faet that the funetion X{t + x) satisfies 

X'{t + x) = f{X{t + x),Ki{X{t + x)),U 2 {t + X-D 2 ),...,Um{t + X-Dm)), 
for all f > 0 and Di < x < D 21 


( 32 ) 
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where we also used the faet that Ui{t + x — Di) = ki (Pi(t + x — Di)) = (X(t + x)), for all x > Di. 
Defining 

P2{D2,t) = PsW, (33) 

and using the faet that p 2 is a funetion of one variable, namely x + t (whieh follows from pT])), one ean 
eonelude that 


P2{t + X-D2) = P2{x,t), Xe[Di,D2]. 


(34) 


Performing the ehange of variables x 
and ( |34l ) we arrive at 


9 + D 2 — t, for all t — D 21 <d<t,m ( [T8] ) and using ( [T^ , ( [28] ), 


P2{e) = Pl{t)+ f f{P2{s),Ki{P2{s)),U2{s),U^{s-D^2),...,Um{s-D^2))ds, 

J t—-D21 

for all t — D 21 < 0 <t. (35) 


Step j: In general, assume that for some j 

Pj (x, t) = Pj{t + X - Dj) 

= X{t + x), X e [-Dj-i,-Dj]. (36) 

We show next that pj+i(x,f) = X{t + x), for all x G [Dj,Dj^i]. We first observe that pj^i{Dj,t) = 
X{t + Dj) = pj{Dj,t). Moreover, the funetion pj+i(x,t) = X{t + x), for all x G [Dj,Dj+i], satisfies the 
following ODE in x 

d^Pj+l{x, t) = f (Pj+l(x, f), Ki (Pj+l(x, f)) , . . . , Kj (Pj+l(x, t)) , Mj+l(x, f), . . . , Mm(x, f)) , 

xG [D,,D,+i], (37) 

sinee the following holds 

X'(t + x) = f {X{t + x),Ki{X(t + x)),..., Kj {X(t + x)), 

Dj+i(f + X — Dj+i),..., Um{t + X — Djn )), for all f > 0 and Dj < x < Dj+i,(38) 

where we also used the faet that Ui{t + x — Dj) = Pi{t + x — Di) = X{t + x), for all x > Dj and i < j. 
By defining 


Pj+i{Dj+i,t) = Pj+i{t), 


( 39 ) 
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once can conclude that 

Pj+i{t + x-Dj+i) = xe[Dj,Dj+i]. (40) 

Performing the change of variables x = 9 + -Dj+i — t, for all t — -Dj+i < 9 < t, we arrive at 

Pj+i{9) = Pj{t)+ j f {Pj+i{s),Ki{Pj+i{s)),...,Kj{Pj{s)),Uj+i{s), 

Jt-Dj+ij 

t/j_i_ 2 (s -Dj_|_ 2 j-i-i); ■ ■ ■) U^i^s diSf for all t Djj^i "Pi 9 "P- t^ (41) 

which completes the proof. 

Note that with this representation we have that 

Uiit) = Ki{j)i{Di,t)), i = (42) 

C. Main Result and its Proof 

Theorem 1: Consider the closed-loop system consisting of the plant ([T3])-([T5]) and the control laws 
( |4^ , (fr7|)-(fT9]). Under Assumptions and there exists a class ICC function (3 such that for all 
initial conditions Xq G M" and Ui^ G C[0,Dj], i = 1,... ,m, which are compatible with the feedback 
laws, the closed-loop system has a unique solution X{t) G (^^[OjOo) and Ui{x,t) G C ([0,Zi)j] x [0, cx))), 
i = 1,..., m, and the following holds 

S(f) < /3 (S(0), f), for all f > 0, (43) 

where 

i=m 

S(f) = |X(f)| + 5^||u,(f)|U. (44) 

i=l 

Corollary 1 (The version ofTheorem^in standard delay notation): Consider the closed-loop system 
consisting of the plant Q and the control laws Under Assumptions and the following 

holds 

i=m / i=m \ 

+ l^*(^)l </S ( l-^(0)l + sup |f/i(6')|,f ) , for alH > 0. (45) 

t-Di<e<t y ^ -Di<e<o J 

The proof of Theorem is based on a series of technical lemmas which are presented next and whose 

proofs are provided in Appendix A. Corollary follows immediately from Theorem by using ([T^. 
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Lemma 1: The backstepping transformations of Wj, i = 1,..., m, defined by 


Wi{x,t) = Ui{x,t) - Ki{pi{x,t)), a: €[0,^1] 


W2{x,t) = U2{x,t) 


K 2 {pi{x,t)), XE[0 ,Di] 
K2{p2{x,t)), XE[Di,D2] 


(46) 

(47) 


t^m{pi{x,t)) , XE[0 ,Di\ 

. . , I f^miP2ix,t)), XE[D^,D2] 

Wm{x,t) = Um{x,t)-< . , (48) 

i,Pm(^X^ ^)) ) X E 

where Pi, i = 1,..., m, are defined in ([T7|)-([T^, together with the eontrol laws ([T7|)-([T9l), transform 
system ([T3|)-([T5|) to the following “target system” 


X{t) = f{X{t),Wi{0,t)+Ki{X{t)),...,Wm{0,t) + Hm{X{t))) 


(49) 


dtWi{x,t) = d^Wi{x,t), xE{0,Di), i = 

Wi{Di,t) = 0, i = 


Lemma 2: The inverse backstepping transformations of (46)-(48) are defined by 


where 


ui{x,t) = wi{x,t) + Ki{'n-i{x,t)), a;e[0,i9i] 


U2{x,t) 


W2{x,t) + 


K2i7Tl{x,t)) , xe[0,T)i] 

f^2{T^2{x,t)) , XE [Di,D2] 


'^m{X) f) 


{ Hm{'Ki{x,t)), xe[0,Di\ 

Km{Tr2{x,t)) , XE [Di,D2] 

. 

iXra^X^ f)) , X E [79m_i, -Dm] 


(50) 

(51) 

(52) 

(53) 


(54) 


7ri(a;,f) 


^it)+ f My, t),wi{y, t) + 1^1 My, t)),---, wM, t) + {-xiiy, t))) dy, 
Jo 


X E [0,Di] 


(55) 


rx 

n2{x,t) = 7ii{Di,t)+ f {7i2{y,t),Ki{n2{y,t)) ,W2{y,t) + K2My,t)), 

J Di 

... ,Wm{y,t) + Umi-xiiyyMdy, xe[Di,D 2 ] (56) 


px 

7 rM,t) = Vrm-l(Dm-l,f) + / f { 7 rm{y,t),Ki{ 7 rm{y,t)) . ,Km-lM{y,t)) , 

Drn — 1 

W^{y,t) + Km{'Xm{y,i)))dy, G [Dm-l,Dm]. 


(57) 
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Lemma 3: There exists a elass JCC function /3i such that the following holds 


2(f) < /3i(H(0),f), for all f > 0, 

(58) 

where 

i=m 

m = ixwi+y^ikiWiioo. 

i=l 

(59) 

Lemma 4: There exist class /Coo functions pi, ..., such that 

lbi(f)l|oo < Pi (2(f)) 

(60) 

lb2(f)||oo < P2(2(f)) 

(61) 

lbm(f)||oo < Pm (2(f)), 

(62) 

where S is defined in (44). 

Lemma 5: There exist class /Coo functions pi, ..., pm such that 

VI 

(63) 

lk2(f)||oo < P2(2(f)) 

(64) 

||^m(^)||oo ^ pm (^(^)) ; 

(65) 

where S is defined in (59). 

Lemma 6: There exist class /Coo functions p, p such that 

2(f) < P(2(f)) 

(66) 

2(f) < p{^it)). 

(67) 

Proof of Theorem^ Combining (67) with (58) we get that S(f) <p 

(/9i (S(0),f)), for all f > 0, and 

hence, with (66) we arrive at ([4^ with (3{s,t) = p{j3i (p(s) ,f)). The proof of existence and uniqueness 

of a solution X{t) G 67^(0, cx)) and ufx^f) G C ([0,i9j] x [0 ,cxd)), i = 

= 1,..., m, is shown as follows. 


Using relations ([20])-([22]) for f = 0, the compatibility of the initial conditions i = 1,... ,m, with 


the feedback laws (42) guarantee that pi{x,0) e C^[Di_i, Di], where Dq = 0. Hence, using relations 
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(p^-([T^, (23H25), and the fact that Ui^ G C[0, A], i = 1,... ,m, it also follows from (46H48) that 
e (^[O, Di\, i = 1,... ,m. The solution to ( [50| ), ( [ST] ) is given by 

D, 




*0 


W, 


( 43 J Wi^{t^x), 0 < X + t < . 

(4.i) = ^o, x + t>D. ■ * = 


( 68 ) 


The uniqueness of this solution follows from the uniqueness of the solution to (50), (511 (see Sections 2.1 
and 2.3 in [dll). Hence, the compatibility of the initial conditions Uig, i = 1,... ,m, with the feedback 


laws (42) guarantee that there exist a unique solution Wi{x,t) G C{[0,Di] x [0, cx))), i = 


From the target system (49) it follows that X{t) G cx)). The fact that 7Ti{x,t) = X{t + x), for 
all t > 0 and X G [i9j_i,i9j], and the inverse backstepping transformations ( |52[ )—( [5^ guarantee that 
Ui{x, t) e C ([0, Di] X [0, CX))), i = 1,... ,m. The proof is completed. 


IV. Stability Analysis Under Predictor Feedback Using Estimates on Solutions 


Theorem 2: Consider the closed-loop system consisting of the plant ([^ and the control laws 
Let Assumptions and hold and assume that the system X = f (X, ki{X),..., is globally 

asymptotically stable. There exists a class XC function /3 such that for all initial conditions Xq G M" 
and Uig G C[—Di, 0], i = 1,... ,m, which are compatible with the feedback laws, there exists a unique 
solution to the closed-loop system with X{t) G U^[0, cx) and A(f), i = 1,... ,m, locally Lipschitz on 
(0, cx), and the following holds 


< /3 (U(0),f), for all f > 0, 


(69) 


where 


U(f) = |x(f)| + V sup \Ui{e)\. (70) 

t-Di< 9 <t 
1=1 — 

Proof: We estimate first |X(f)|, for all f > 0. Since the system X = / (X,a;i,... ^Um) is forward 
complete using Lemma 3.5 from [l25l and the fact that /(O, 0,..., 0) = 0 (which allows us to set R = 0), 
we get that 


|X(f)|<A l^(0)| + V sup |A(0)| , forall0<f<A, (71) 

V ^^-Di<e<o j 

for some class /Coo function A- Using the fact that Ui{t) = ki (Pi(f)), for all f > 0, with Pi{t) = 
X {t + Di) and the fact that system X = / (X, ki{X),u: 2 , • • •, ccm) is forward complete with respect to UJ 2 
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we get by applying Lemma 3.5 from [1251 that \X{t)\ < ^>2 + EL™ ^^PDi-Di<e<D 2 -D, \Ui{0)\y 

for all Di ^ t ^ D 2 t for some cl3.ss KLoo function ' 2 /^ 2 ’ Hence, since Dj ^ Di, Vj < i, with ( [tTI we get 

( i=m \ 

|X(0)| + V sup \Ui{e)\ , for all < t < D^, (72) 

7^ -Di<e<o J 

where the elass /Coo function V’ 2 (s) is defined as -02 (s) = L 2 (f/’i('S) + s). By repeatedly applying Lemma 
3.5 from [|251 we get under Assumptionthat there exist class /Coo functions ipj, j = 3,... ,m such that 
|X(/)| < Li [\X{Dj_i)\ + ELr \Uii0)\), for all Dj_i < t < Dj, and hence, 

( i=m \ 

|X(0)| + V sup |f/L^)l , for all (73) 

-D,<e<o J 

where the class /Coo functions 'ifjj, j = 3,..., m are defined as 'ipj{s) = rjjj ('0i-i(s) + s), j = 3,..., m, 


and where we also used the fact that Dj < Di, Vj < i. Combining (71 )-(73) we arrive at 


'3 

i=m 


|X(/)|<^j |X(0)| + V sup \Ui{e)\], for all 0 < / < 

' —' -Di<e<o ' 


(74) 


2=1 


where 'ij(s) = Using the fact that Umit) = Km{Pmit)), for all / > 0, with Pmit) = 

X {t + Dm) and the fact that system X = f (X, ki{X), . .., Km{X)) is globally asymptotically stable we 
get that |X(/)| < $i{\X{Dm)\,t- Dm), for all t > Dm, for some class ICC function /3i. Hence, using 
( |74| ) we get that 

|X(/)| </32(U(0),/), for all/>0, (75) 

where the class JCC function ^2 is given by (32{s,t) = L (-0(5), max {0, / — Dm}) + 

We estimate next sup^_JJ_^^g^^\Ul{6)\. Since Ki is locally Lipschitz and Ki(0) = 0, there exists a 
class /Coo function di such that |ki(X)| < di (|3C|), for all X G M”. Since for all / > 0 it holds that 


Ui{t) = Ki {X{t + Di)) using (75) one can conclude that 


sup \Ui 
t-Di<e<t 


< sup If/L^)! + di (/32 (U(0),0)) , 0<t<Di, 

-Di<e<o ^ ^ 


and hence. 


sup |Ui(6*)| < d 2 (U(0)), 0 < / < D\, 

t-Di<0<t 


(76) 


(77) 


where the function d 2 (s) = di (L 2 (s, 0)) + s belongs to class /Coo- Using identical arguments we also 
get that 


t-Di<e<t 


sup < di (Ls (U(O),/:)) , t>Di 


(78) 
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Combining (|77|) with (78) we arrive at 


sup |?7i(6')| < /^s (fi(0), t), for alH > 0, (79) 

t-Di<e<t 

where the elass /C£ funetionis defined by/33(s, f) = di (^^2 (f2(0),max{0,t — 

Using the faets that K 2 is loeally Lipsehitz and that ^2(0) = 0, whieh allows one to eonelude that there 
exists a elass /Coo funetion 03 sueh that |k 2(-^)| < “3 (|X|), for all X G M”, with similar arguments we 
get that 


sup |U2(6*)| </34 (U(0), f), for all f > 0, 

t-D2<e<t 


(80) 


where the elass /C£ funetion/34 is defined by/34(s, f) = ds (^^2 (U(0),max{0,t — i92})j+Q;4(s)e“™'^^'f°’*“^2^, 
where the funetion 04(5) = da (^^2 (s, 0) j + s belongs to elass /Coo- With the same arguments one ean 
eonelude that there exist elass ICC funetions /3j+2, j = 3,..., m, sueh that 


sup \Uj{9)\ < ^j +2 (^(0), t), j = 3,...,m, for all f > 0. 


(81) 


Combining estimates (75), (79), (80), and (81) we get (69) with (3{s,t) = (32{s,t) + A+2('5,/)• 

From ([^ using the faet that Ui^ G (C[—i9j,0], for all i = the Lipsehitzness of the veetor 

field / guarantees that X{t) G C^[0, Di). The faet that Ui{t — Di) = ki {X{t)), for all t > Di, and the 
Lipsehitzness of ki guarantee that X{t) G C^{Di, D 2 ). Sinee Ui^ is eompatible with the first feedbaek law 
one ean eonelude that X{t) G (7^ [0,792). Analogously, sinee for f > 79i the state X evolves aeeording 
to X{t) = f {X{t), Ki {X{t)) ,U 2 {t — D 2 ),... ,Um{t — Dm)), the faet that U 20 is eompatible with the 
seeond feedbaek law and the faet that U 2 {t — D 2 ) = ^2 iX{t)), for all t > D 2 , where k ,2 is loeally 
Lipsehitz, guarantee that X{t) G C^[0, D^). Continuing this proeedure it is shown that X{t) G C^[0, cx)). 
Sinee Ui{t) = Ki {X{t + Di)), i = 1,... ,m, the Lipshitzness of Ki, i = 1,... ,m, guarantees that Ui{t), 
7 = 1,..., m, are Lipsehitz on (0, 00 ). ■ 

V. Stabilization of the Nonholonomic Unicycle Subject to Distinct Input Delays 


We eonsider the following system 



= U2(f- 792 ) COS ( 963 ( 7 )) 

(82) 

Mt) 

= f/2(7-I92)sin(X3(7)) 

(83) 

L(«) 

= f/i(7-I9i), 

(84) 
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which describes the dynamics of the unicycle, where (Xi,X2) is the position of the robot, X3 is heading, 
f/2 is speed, and Ui is the turning rate. We consider the following time-varying nominal control law 
designed in 


Ui{t) = — cos (f) — M(f)Q(f) (1-f cos^ (f)) — X3(f) 
U2it) = -M{t) + Q(t) (sinit) - cos(t)) + Q(t)Ui(t), 


( 85 ) 

( 86 ) 


where 


M(f) = Xi(f)cos(X3(f))+X2(f)sin(X3(f)) 
Q{t) = Xi(t)sin(X3(t)) - X2(t)cos(X3(t)), 


( 87 ) 

( 88 ) 


which achieves global asymptotic stabilization when Bx = £>2 = 0 . We employ next our predictor- 
feedback design when 0 < £>1 < £>2- We first verify that Assumptions and are verified for system 


( 82 H 84 ) under the control laws ( 85 H 88 ). We first note that system X = / (X, a;i,a;2), where 


/ (-^5 ^1, ^2) — 


022 COS (X3) 
0J2 sin (X3) 
021 


( 89 ) 


is forward complete with respect to (a;i,a;2), and hence. Assumption is satisfied. One can see this by 
defining 


i=3 


) = 5 E 


(90) 


i=l 


and readily verifying by employing Young’s inequality that 

dR(X) 


dX 


- f (X, 021,022) = XicC2Cos(X3)-i-X2a;2sin(X3)-f X3CC1 
< R(X) + -{021 + 02I) . 


( 91 ) 


Moreover, under ( 85 ) the dynamics of the nominal, delay-free system are described by X = / (X, ki (t, X ), 022), 
where 


f(X,nx (t,X),022) = 


022 cos (X3) 
022 sin (X3) 


_ -X3 - M (Xf cos (t) -M(X)Q (X) (1 -f cos 2 (t)) 
Therefore, by defining S'2(X3) = |X| we have that 

^^ 2 (^ 3 ) 


( 92 ) 


0X3 


(-X3 - M(Xf cos (t) - M(X)Q(X) (1 + cos' (t))) < 4 (Xf + X^) , 


( 93 ) 
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where we employed Young’s inequality and the faet that = Xf + whieh follows from (87) 


and (g. With ^i(Xi, X 2 ) = ^ {Xf + X|) we get that 

a^i(Xi,X2 


dXi 


, dS^(X,,X2) . , 1 4 

U2 cos (X3) H-——- U2 sm (X3) < 3 oi + -UJ2- 

UA.2 ^ 


(94) 


Henee, defining S = Si + S 2 for system X = f (X, ki (f, X), U 2 ) with (92) it follows that 

—^^^f{.X,Ki{t,X),U2) < IQS'! + -072 

< IQS' + -a; 2 - 


(95) 


Sinee the differential inequality for S (X(f)) along the trajeetories of the system is linear in S (X(f)) 


it follows that system X = / (X, ki (f, X), 012 ) with (92) is forward eomplete with respeet to 002 - 
Therefore, Assumption hold^ The elosed-loop system in the ease Di = 792 = 0 is uniformly globally 
asymptotieally stable as it is proved in [l49ll . and henee, one ean apply. Theorem 
The predietor-feedbaek eontrol law is given by 


Ui{t) = -Mi{tycos{t + Di)-Mi{t)Qi{t){l + cos^t + Di))-Pi^^{t) (96) 

U 2 {t) = -M 2 {t) + Q 2 {t) {sin{t + D 2 ) - cos{t + D 2 )) + Q 2 {t) 

X {-M 2 {tf cos (t + D 2 ) - M 2 {t)Q 2 {t) (1 + cos^ (t + 792)) - P2X3 if)) , ( 97 ) 

where 


Mif) = it) cos {Pi^^ (t)) + Pi^^ (t) sin {Pi^^ (f)) , i = 1,2 


(98) 


Qif) = Pi^^ f) sin [Pi^^ {t)) - Pi^^ {t) cos {Pi^^ {t)) , 7 = 1,2. 


(99) 


* It can be shown that M{t) = Uf{t — D 2 ) — Q(t)Ui{t — Di), Q{t) = M{t)Ui{t — Di), and X 2 ,{t) = C/i(i — Di). Therefore, when 

' U 2 {Y) - Y 2 UJ 2 ' 

to hold, the system Y = f {Y, ki (t, Y) ,ui 2 ) with / (Y, ki {t, Y ), 11 ) 2 ) = 


Di > D 2 in order for Assumption 


Y 1 W 2 

UJ2 


where 


U 2 {Y) = -Y 1 +Y 2 (sin(t) - cos(tj)-YiX 2 cos{t)-YiYi (l + cos^ (f)) — ^ 2 ^ 3 , has to be forward complete with respect to 022 - However, 
this might not be the case. Consider, for example, the case in which D 2 = 0, cu 2 = 0, ^ 2 ( 0 ) = 1, ^ 3 ( 0 ) = 0, and hence, Y 2 {t) — 1, for all 
t > 0. The state Yi satisfies Yi = — Yi + sin(f) — cos(f) — Y^ cos (f) — Yi (l + cos^ (t)), for all f > 0. We show next that Yi escapes to 
infinity before t ~ j- Choose Yi(0) = —fi < 0, where /3 is sufficiently large such that Yi(f) < 0 for all f As long as Yi(f) < 0 and 
^ < f it holds that Yi < —3Yi — Yi^|. Hence, using the comparison principle we get that Yi(f) < — ■ Choosing, for example. 


/3 = 12 we have that Yi(f) < — ■ Hence, Yi(f) < 0 for all t < 

(which also implies that Xl —>■ 00 ). 


l°g(2) 

3 


0.23 < f « 0.79. Moreover, |Yi| 


00 before t = 
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The predictors are given by 

= x,{t)+[ U2{e-D2i)cos{Pi^^{e))de 

Jt-Dx 

(t) = X 2 {t) + [ U 2 {e - D 21 ) sin (Pi {6)) d9 


' t — D\ 

rt 


Pi (t) = X3(t)+ / U,{9)d9 

^ Jt-Di 


and 


rt 


' 4—D 21 
rt 


P2xS^) = -Plxi(^) + 
-P2X2 (^) = P^X 2 (^) + 
-P2X3 (^) = -Plxg (^) + 


't—D21 


x.r, - -lx.r, . I £/ 2 ( 9 ) cos (P 2 x, («))<*« 

t—D21 
rt 

U2{9) sin {P2x^{9))d9 
k,{9 + P 2 , P 2 ,,, (0), P2x 2 {0). P2x, md9, 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 


where for all t — P 21 ^ 9 <t 

K,{9 + D2,P2xSd),P2x2id),P2x,m = -M2{9 YcOs{9 + D2)-M2{9)Q2{9) 

X (1 + cos" {9 + D 2 )) - P2x,{0)- 


(106) 


Note that the Pi-time units ahead predictors, namely Pi^^ , i = 1,2,3, given by (100H102) can be 


computed explicitly in terms of the history of Ui{9) on the interval 9 G [t — Di,t], the history of 1 / 2 ( 9 ) 
on the interval 9 G [t — D 2 , t — P 21 ], and the current states Xi(t), i = 1, 2, 3, as 

Pi^^(t) = Xi(t)+ [ U2(9-D2i)cos(x3(t)+ [ Ui(s)ds]d9 (107) 

Jt-Di \ Jt-Di J 

Pi^(t) = X 2 (t)+ [ U2(9-D2i)sm(x3(t)+ [ Ui(s)ds]d9 (108) 

Jt-Di \ Jt-Di J 

Pi (f) = X3(f)+ f U,(9)d9. (109) 

® Jt-D. 


This is not possible for the P 2 -tlnie units ahead predictors, namely P 2 x , i = 1, 2, 3, given by (103)-(1061, 


since in this case P 2 x defined in (105) can not be solved explicitly in terms of the current states of the 


plant and the history of the actuator states. 

VI. Application to Linear Systems 

We specialize our control design to the case of linear systems in which case the control laws are given 
explicitly. We consider the following system 

i=m 

X(t) = AX(t) + Y,W(t-D,), 

i=l 


( 110 ) 
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which can be written as 


X{t) 

dtUi{x,t) 
Ui {Di,t) 


AX{t) + y^^biUi{0,t) 

i=l 

d:,;Ui{x,t), xe(0,A), i = 


( 111 ) 

( 112 ) 

(113) 


where ^4 is an n x n matrix and bi, i = 1 ,... ,m are n-dimensional vectors. System ( |111[ )-( [TT3 1 is the 
linear version of system ([T3|)-([T5|). We assume linear nominal feedback control laws, that is, in the delay- 
free case we have Ui(t) = kfX(t), and hence. Assumptionis satisfied with Ki{X) = kfX,i = 1,..., m, 
under the assumption that the pair (A, [ bi ... 6^ ]) is stabilizable. Note that Assumptions andhold 
for the case of linear systems under linear nominal feedback controllers. We first re-write the predictor 
states in their PDE representation, namely (20)-(251, for the case of linear systems given by (111 )-( |l 13| ) 
as 


d^pi{x,t) = Api{x,t)+ '^biUi{x,t), xe[ 0 ,Di] 

i=l 

i=m 

dxP 2 ix,t) = {A + bikl) P 2 {x,t)+ '^biUi{x,t), x e [Di, 02 ] 


(114) 

(115) 


2=2 


i=m—l 


dxPm (^7 1 ^ ^ ^ bikj^ I “1“ bjYiUm t), X G Dm\ 7 


(116) 


2=1 


with initial conditions given by (23)-([25]). Solving explicitly the linear ODEs in x (114Hi 16) and using 


the boundary conditions ([23])-( 25) we get that 


Pi{xA) = e^'^X{t)+ y'^'^biUi{y,t)dy, x e [0, E>i] (117) 

do i=i 

2—772 

P2{xA) = e'^d^~^^'^pi{Di,t) + xe[A,A] (118) 


i=2 


Pm{x,t) = b^aUmiy, t)dy, 

'X Dm — 1 


X G [-D772—17-^m] 7 


(119) 
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where we used the notation 


An = A 


]=i 


The eontrol laws are given by 


i=i 


Ui{t) = Kjpi{Di,t), i = 


( 120 ) 

( 121 ) 


( 122 ) 


Theorem 3: Consider the elosed-loop system eonsisting of the plant (|1 1 1|)-([TT3) and the eontrol laws 


(122), (117)-(119). Let the pair (A, [ bi ... 6m]) be stabilizable. There exist positive eonstants A and 
p sueh that for all initial eonditions Xq G and Ui^ G 77^(0, Di), i = 1,..., m, whieh are eompatible 
with the feedbaek laws, the elosed-loop system has a unique solution (X(f), mi(-, f),..., Mm(-A)) G 
C ([0, cx)); X i7^(0, Di) x ... x if^(0, Dm)) ([0, ex); x L^(0, Di) x ... x L^(0, Dm)) and the 
following holds 


T{t) < nT{0)e for all f > 0, 


(123) 


where 


i=m „£) 


T(t)^\X{t)\ + Y. dx . (124) 

i=i do 

Note that with definitions Pi{x,t) = Pi{t + x — Di), for all x G [A-i, A], with Dq = 0 (see Seetion 
III-B) the predietors (117 )-( TT9] ) ean be written as 

i=i 

biUi{s — Dii)ds, t — Di < 9 <t 


i=m 

P^{e) = 

Jt-Di 


P^{e) = + / e 

J t — D2l 


2=1 

Q i=m 

Aiie-s) 


'^biUi{s - Di2)ds, t-D2i<9<t 


(125) 

(126) 


i=2 


Pmie) = / e^-^^^-%mUmis)ds, t - Dmm-1 < 9 < 1,(127) 

^ Dmm — 1 


where the matriees Aj, i = — 1, are given in (120), (121). With this notation the eontrol laws 


(122) are written as 


Ui(t) = KjPi(t), i = 


( 128 ) 
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We have the following corollary as an immediate consequence of Theorem and relation ( [T^ . 

Corollary 2 (The version of Theorem^in standard delay notation): Consider the closed-loop system 


consisting of the plant (110) and the control laws (128), (125)-(127). Under the assumption that the pair 
{A, [hi ... bm ]) is stabilizable the following holds 

i=m / i=m „o \ 

/ Uf&fde < /i |X(0)| + V / Ui(efde e-^\ for all t > 0. (129) 

i=i \ i=i J-D. ) 

Proof: We prove Theorem by showing that the control laws (122), (117)-(119) are identical to the 
ones introduced in [f53]l and [|54ll . and hence. Theorem is proved by following the proof of Theorem 1 
in dm. Equivalently we show that the backstepping transformations ([^ specialized to the linear case are 
identical to the backstepping transformations introduced in dSl- Toward that end define 

(130) 


/ \ \ X, X < Di 

- I x>Di ’ « - 

Pi{xA), xe[0,Di] 

, ,, , P2ixA), xe [Di,D2] . 

9i{x,t) = < . , 1 = 1, ...,m. 

Pi{x,t), xe[Di_i,Di] 


(131) 


Then, using (117)-(119) it follows that 


e ^'-'^'^gi{x,t) = Vi_i{x,t), xe[ 0 ,A], i = l,...,Tn, 


(132) 


where for all a: G [0, Dm] 

Vi{x,f) = f) + / hjUj{y,t)dy 

j=i+i 

^01 (x) j=m 

Vo{x,t) = X{t)+ e~^^^bjUj{y,t)dy. 

Jo j^i 

We show this by induction. For all x G [0, Df it holds that 


(133) 

(134) 


e (x,f) = e ^'"pi(x,f) 


—Ax^ 


= m + 

= Vo{x), 




^-Ay 


^bjUj(y,t)dy 


i=i 


(135) 
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where we used ( |117[ ) and ( |130| ). Assume now that ( |132[ ) holds for some i. It follows that 


e = e { 


Ai‘, 


= e 


' Pl{x,t), 

X G 

[o,A] 

P2{x,t), 

j 

X G 

[A, A] 

\ 

X G 

[A-i, A] 

, Pi+l{x,t), 

X G 

[A, A+i] 


x,t), 

X G [0, A] 

1 Pi+iix,t), 


X G [A, D 


j+ij 


for all X e [0, A+i]- Using (117Hi 191 and (130) we get from (136) that 

_^,3, f f), xe[0,A 


e '^gi+i{x,t) = 


e + /^, e EL™ i bjUj{y, t)dy, x e [A, A+i] 


Using ( |131| ) and ( |132[ ) it follows that 

and henee, 


(136) 


(137) 


(138) 


Pi{Di,t) = e^‘-i"Hi_i(A,t) 


(139) 


Combining (137) with (139) we arrive at 


.-Ax, / _ r e _ xe[0,A] 

’ \ e + /oA [A,A+h ' 


We then observe from (130) that (j)i{x) = (j)i+i{x) = x, for all x < Di, and henee, 

/.0,+i(x) i=™ 


e ^ bjUj{y,t)dy = 0, x e [0, A]- 

j=i+i 


(141) 


Moreover, using ( |130| ) we get that 

g-AiDigA-iA ^ g-A0i(a^)gA-i'^'U)^ xe[A,A+i]- (142) 

The proof is eompleted if we show that 

Vi-i{Di,t) = Vi-i{x,t), for all a; e [A,-Dm] and i = 1,... ,m. (143) 

We prove this by induetion and by using definitions (133), (134). Using (|130) for z = 1 we have that 


VoiDi,t) =X{t)+ e 

Jo 


Di 

-Ay 


^bjUj{y,t)dy 

i=i 


/•(^i(x) i— 

= X{t)+ e~^'^ bjUj{y,t)dy, for all x > Di 

Jo 

= no(a;,f), for all x > Di. 


(144) 
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Assume now that (143) holds for some i. Then using (130) and (133) we have that 

/•A+i 


Vi{Di+xA = e + / e 

Jd, 


-A; 


bjUj{y,t)dy 

j=i+l 

p4ii+i{x) 

j=i+i 

for all X > A+i- 


(145) 


Therefore, using ( |143| ) we get that 

P<j>i+l{x) A™ 

Vi(A+i,t) = ^_ / g-Ajy bjUj{y,t)dy 

j=i+i 

= Vi{x,t), for all a: > A+1, (146) 


whieh eompletes the proof that (132) holds. Therefore, using (131) the baekstepping transformations 


(46)-(p8l) ean be written in the linear ease as 


Wi{x,t) 

= Ui{x,t) — Kje'^°^Vo{x,t), 

X G [0,T>i] 

(147) 

W2{x,t) 

= U2{x,t) — K2e^^^Vi{x,t), 

X G [0, D 2 ] 

(148) 

Wm{x,t) 

= Um{x, t) - 

(x,f), xG[0,i7m], 

(149) 


whieh are identieal to the baekstepping transformations (67), (68) from 


As a speeial ease we provide explieitly the first two baekstepping transformations. Using (133), (134) 


the baekstepping transformations (147), (148) take the form 

/ PX 

Wi{x, t) = Ui{x, f) - fcf ^ kuiiy, t)dy j , a; G [0, Di] 


(150) 


i=l 


kj (e^^X{t) + £ y'^ biUi{y, t)dydyj , x G [0, Di] 

W2{x,t) = U2{x,t) — < /cj biUi{y,t)dy ,(151) 

+ /di EC™ biUi{y, t)dy^ , x G [T>i, A] 

whieh are identieal to relations (15), (77), and (78) from [[54l . The predietor feedbaek eontrol laws are 




ui{Di,t) = fcf (^e^^^X(f) + y y'> {biui{y,t)+ b2U2{y,t))dy 

U2{D2,t) = fcj y e^^^^~y^biUi{y,t)dy 

_(_gAiX>2i f e'^^^^~y^b2U2{y,t)dy + f e'^^’'^^~y'>b2U2{y,t)dy\ , 

Jo Jdi J 


(152) 


( 153 ) 
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and in terms of the delayed states Ui{t + x — Di) = Ui{x, t) by 

rt rt—D 2 i 


/ ft ft-U2\ \ 

Ui{t) = kUe^^^x{t)+ e^^*-^'>hUi{e)de+ e^^^-^-^^^^b2U2{e)de) 

V Jt—Di Jt — D 2 / 

U2{t) = + [ e^^^-^^biUi{e)de 

V Jt-Di 

ft—D2\ ft \ 

+ / e^^^-^-^^^\U2{0)de + / e^^^^-%2U2{0)de . 

Jt—D2 Jt — D21 / 


(154) 


(155) 


Relations (154), (155) ean be written as 


U,{t) = fcfPi(t) 


U2{t) = ki + [ e^^^^-%2U2{e)de] , 

V Jt-D21 J 


(156) 

(157) 


whieh are the eontrol laws (46), (47) from [|54l . Sinee the explicit expression of the general i-th control 
law for the general case of m inputs is identical to the one obtained in ll54l (relation (43)), for clarity of 
exposition we provide it again here 


rt 


u^it) = kj <h(A,o)x(f) + ^ (A, e-t + Dj)b,u,{e)de 


j=i 


3=m 


E 


^{Di,e-t + Dj)bjUj{e)de , 


(158) 




where 


<h(a;,|/) = Di)^Ai.j_(Di Di_i) ^ ^ ^ ^Aj{Dj+i y) ^ Di < X < Di+i and Dj < y < (j)j+i{x),(l59) 


for any i,j G { 0 , 1 ,..., m — 1 } satisfying 0 < i < j, and 


^{x,y) = for all < x < y < Pj+i, 


(160) 


for i G {0,1,..., m — 1}, where (j)i is defined in (130) and we adopt the notation Pq = 0 


VII. Simulations 

We consider the stabilization problem of the nonholonomic unicycle subject to distinct input delays 
from Section [V| whose dynamics are described by (82)-(84). We choose Di = 0.5 and D 2 = 1. The 
initial conditions for the plant are chosen as 7fi(0) = ^ 2 ( 0 ) = ^ 3 ( 0 ) = 0.5 and for the actuator states as 
Ui{d) = 0, for all —Pi < 6 ^ < 0, and U2{0) = 0, for all —P 2 < 6 ^ < 0. In Fig. we show the response of 
the closed-loop system under the predictor-feedback controller (96)-( 106) and in Fig. the corresponding 
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Fig. 3. The closed-loop response of system l|82|l-l|84^ with Di = 0.5 and D 2 = 1 under the predictor-feedback controller j96|l- i fToel l. The 
initial conditions of the plant are Xi(0) = X2(0) = X3(0) = 0.5 and for the actuator states are Ui{8) = 0, for all —0.5 < 0 < 0, and 
U 2 {e) = 0 , for all -1 < < 0. 


control efforts. The predietor-feedbaek eontroller asymptotieally stabilizes the nonholonomie unieyele. In 
partieular, at f = 1 both delays are eompensated and the system behaves as in the nominal, delay-free 
ease. In eontrast, at is is shown in Fig. the elosed-loop system under the nominal, uneompensated 
eontroller is unstable. 

VIII. Conclusions 

We presented a predietor-feedbaek eontrol design methodology for the stabilization of multi-input 
nonlinear systems with distinet input delays. We proved global asymptotie stability of the elosed-loop 
system using Lyapunov arguments and arguments based on estimates of elosed-loop solutions. We also 
dealt with linear systems as a speeial ease of our methodology. We applied our approaeh to the stabilization 
of the nonholonomie unieyele with delayed inputs. 

In eontrast to the single-input ease, for whieh the predietor-feedbaek eontroller is available explieitly 
for some elasses of nonlinear systems (see, for instanee, [|27l ) with a speeifie open-loop strueture, in 
the multi-input ease the elass of nonlinear systems for whieh the predietor-feedbaek eontrol law ean be 
obtained explieitly seems to be more restrietive (although it can be obtained explieitly in some trivial 
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Fig. 4. The control effor ts l|96| l and \91\ of the closed-loop response of system l|82[(-l|84[( with Di = 0.5 and D 2 = 1 under the predictor- 
feedback controller i96i-il06i. The initial conditions of the plant are Xi(0) = X2(0) = ^3(0) = 0.5 and for the actuator states are 
Ui{e) = 0, for all -0.5 < 0 < 0, and 1/2(61) = 0, for all -1 < 61 < 0. 



Fig. 5. The closed-loop response of system (|82^-(|84^ with Di = 0.5 and D 2 ~ 1 under the nominal controller l|85[l-j88[l. The initial 
conditions of the plant are Xi(0) = 2^2(0) = X 3 J 0 ) — 0.5 and for the actuator states are U\{6) = 0, for all —0.5 < 0 < 0, and U 2 (d) = 0, 
for all -1 < 6» < 0. 
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cases, such as, for example, the case of linear systems). This is attributed to the faet that the formulae of 
the predictors in the multi-input case depends not only on the open-loop structure of the system but also 
on the form of the feedbaek functions. 

As a next step, it is of interest to study the problem of stabilization of multi-input nonlinear systems 
with actuator dynamics governed by wave or diffusion PDEs with different wave propagation speeds or 
diffusion coefficients, respeetively, in eaeh individual input ehannel. The starting point for sueh a study 
is ttH. 


Appendix A 


Proof of Lemma 


Setting a; = 0 into (46H48), and ( [Tt] ) we get (49). Using (26) one ean eonclude that dtpi{x,t) = 
dxPi{x,t), for all x G [0,Ui]. Henee, using ([T^ and (46) we get that 


dtWi{x,t) - dxWi{x,t) = dtUi{x,t) - dxUi{x,t) -f _ d^pi{x,t)) 


dp 


= 0 . 


(A.l) 


Similarly, using the fact that pfx^t) = X{t -f x), for all x G [A-i, A] and i = 2,...,m, we get (50). 
Setting X = Di into (|46l)-(|48]) and using ( [T?] ), ( [42] ) we arrive at ( [5T] ). 


Proof of Lemma 

We prove this lemma by showing that pi{x,t) = 7ri(x,f), for all x G [0,i9i], and Pi{x,t) = Trfx^t), 
for all X G Di], i = 2,... ,m. Equivalently we show that 'xfx, f) = X{t + x), i = 1,..., m and 

use the faet that Pi{x,f) = X{t + x), i = 1,... ,m. We first observe that 7ri(0,f) = X{t), and henee, tti 
satisfies the following initial value problem for all a: G [0, Ui] 


dx'ni{x,t) = f {'Ki{x,t),Wi{x,t) + Ki{ji{x,t)) , . . . ,Wm{.X,f) + rirn{xi{x,t))) (A.2) 

7ri(0,f) = X{t). (A.3) 


The solution to (A.2), (A.3) is Tii{x,t) = X(t + x). This solution satisfies the boundary eondition (A.3). 


It also satisfies the ODE (A.2) since by (49), (50) it follows that the following holds 


X'{t + x) = f {X{t + x),Wi{t + x-Df + Ki{X{t + x)), 

..., Wm{t + X — Dm) + Km {X{t + x))), for all f > 0 and 0 < x < Di, (A.4) 
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where the solutions to ( |^ , ( [ST] ) are given by Wi{x, t) = W{t-]-x—Di), for all x G [0, Di] and i = 1 ,..., m, 
where Wi, i = 1,... ,m, satisfy Wi{t) = 0, for all f > 0. Assume now that Trj{x, t) = X{t + x), for all 
X e [Dj_i,Dj], and some j. Then we elaim that vr^+i = X{t + x), for all x E [Dj,Dj+i\. The funetion 
TTj+i satisfies the following initial value problem 


d^7lj+l{x,t) 


f {7ij+i{x, t), {nj+i{x, t)) {nj+i{x, t)) , Wj+i{x, t) 


'^j+l y 


+ Kj+i{'Kj+i{x,t)) , . . . + K^{7lj+l{x,t))) , XE [Dj,Dj+i] (A.5) 

TTj{Dj,t). (A. 6 ) 


Sinee TTj{x, t) = X{t + x), for all x E [-Dj-i, Dj], it follows that 7rj+i(x, t) = X{t + x), x E [Dj, -Dj+i], 

satisfies the boundary eondition (A. 6 ). Using ( [50] ), ( [ST] ) we get that Wi{x,t) = Wi{t + x — Di), for all 
X E [0, Di], where Wi{t) = 0, for all f > 0. Sinee from ([4^ it holds that 


X'{t + x) = f{X{t + x),Wi{t + x-Di) + K,{X{t + x)), 

..., Wmit + X — Dm) + Km {X{t + x))) , for all f > 0 and a: > 0, (A.7) 


one ean eonelude that the following hods for all f > 0 and x E [Dj, Dj+i] 

X'{t + x) = f {X{t + x), Ki{X{t + x)),..., Kj {X{t + x)) ,Wj+i{t + X - Dj+i) 

+ Kj +1 {X[t + x)) , . . . , Wm{t + X — Dm) + Km {X {t + x))) , (A.8) 


and henee, tTj+i = X{t + x), for all x E [Dj, Dj+i]. Sinee this holds for an arbitrary j one ean eonelude 
that TTj = X{t + x), for all x E [Dj_i,Dj] and j = 1,..., m, with Do = 0, whieh eompletes the proof. 


Proof of Lemma 

From Assumption it follows from [153 that there exist a smooth funetion S (Z) : —)■ M_,_ and 

elass /Coo funetions ai, 02 , 03 , and 04 , sueh that 


dS (X) 
dX 


«i(|X|) < 5(X) <«2(|X|) 

/ (-A, Ki (X) + (Ui, . . . , Km (X) + bJm) < “^3 (1^1) + j j ) 


(A.9) 
(A. 10) 


, 2 = 1 


for all X G M" and cui,... ,0Jm ^ IR- With similar ealeulations as in 11371 (Theorem 5) we get that 

d\\wi{t)\\c, 


dt 




(A.ll) 
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along the solutions of ( [50{ ), ( [5T] ). Consider the Lyapunov funetional 

2 /-E-ir IKWIIc,oo 

v(t) = S(x(t)) + t / 


(A.12) 


Using Lemma (relation (|^) and relations ( |A.10| ), ( |A.11[ ) we get along the solutions of (|^-([5T|) that 

( i=m \ / i=m 

^|wi(0,f)|j -2^4 ( ^ ||wi(f)||c,oo 1 • (A.13) 

Using the faet that |tej(0,f)| < ||tej(f)||c,oo, i = 1,..., m, we get that 

/ i=m 

V{t) < -«3 (|X(f)|) - «4 \\w^mc,^ 1 . (A.14) 


, 2 = 1 


It follows with the help of (A.9) that there exists a elass K, funetion 71 sueh that 


rw <-71 (i^w). 


(A.15) 


and henee, with the eomparison prineiple (see, for example, ll33ll l one ean eonelude that there exists a 
elass )CC funetion /?2 sueh that 


V{t)<(32{V{0),t). 


(A.16) 


Using (A.9), the faet that ||tej(f)||oo < \\wi 


< e^^^\\wi{t)\\oo, i = 1 ,... ,m, and the properties of 


elass 1C funetions we get estimate (58). 


Proof of Lemma 


We prove the lemma by induetion. For elarity we present two initial steps. We prove first bound (60). 
Using the faet that pi satisfies (20) we get under Assumption that there exists a smooth funetion 
R\MP ^ M+ and elass /Coo funetions 05 , ag, and aj sueh that 


dR{X) 

dX 


« 5 (|X|) < R{X)<a^{\X\) 

/ i=m 

f{X,uj,,...,uJ < R{X) + arlj2 


UJi 


(A.17) 

(A.18) 


. 2 = 1 


for all X G M" and ce* G M, i = 1,..., m. Therefore, 
dR{pi{x,t)) dR{pi{x,t)) 


dx 


dpi 


f iPi{xR),Ui{xR ),... ,M^(a;,f)) 


< 7((pi(a;,f))-f Q!7 ( ^ |Mi(a;,f)| I , xe[0,Di]. 


(A.19) 


, 2=1 
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,Di 


Hence, using ( [23] ) we get that 

R{pi{x,t)) < R{X{t)) + J a-i ^^^^\ui{x,t)\ ^ dx, (A.20) 

and hence, 

( i=m 

^ll^iWIloo 1 . (A.21) 

With the help of ( jA.lVj ) and the properties of class /C functions we get estimate ( [60| ). We prove next ( [M] ). 
Under Assumption one can conclude that there exists a smooth function : M” —)■ M+ and class /Coo 
functions ag, ag, and aio such that 


asd^l) < R^{X) < ag{\X\) 
—^ gi{.X,uj2, ■ ■ ■ ,ujm) < Ri{,X) + aiQ y 


Ui 


(A.22) 

(A.23) 


A=2 


for all X G M” and cuj G M, / = 2,..., m, where gi is defined in Assumption Using the fact that p 2 
satisfies ( pTj ) we get from ( |A.23| ) that 

dRi {p2{x, t)) dRi {p2{x, t)) 


dx 


dp 


-9l iP2{x, f), U2{xR),..., Umix, t)) 


< (P2(a;,t)) + aio ( 5]] I , a: G [U>i, £> 2 ]- 


(A.24) 


J=2 


Thus, using (24) we get that 


-02 


Ri{P 2 {xR)) < e^^Ri{pi{DiR)) + aip \ui{x,t)\ dx, 


(A.25) 


, i =2 


and hence, from (A.22) that 


Ri{p2{x,t)) < e^^ag{\\pi(t)\\oo) + D2e^^aio i'^\\ui(t)\\, 


(A.26) 


A=2 


Using (A.22) and (60) we get from (A.26) estimate ([M]). Assume now that for some j it holds that 


IbiWIloo < Pj (s(f)) 


(A.27) 


Under Assumption there exist a smooth function Rj : M” —)■ M+ and class /Coo functions Q! 8 + 3 (i-i)> 
Q^ 9 + 3 (i-i)> ^nd aio+ 3 (j-i) such that 


dRA^) 

dX 


tt8+3(j-i)(|2C|) < Rj{X) < a9+3Q-_i)(|X|) 

( i=m 

E \x!i 

i=j+l 


(A.28) 

(A.29) 
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for all X G and a;* G M, i = j + 1,... ,m. Using the fact that pj+i satisfies (37) and definition 


gj (X, ccj+i,..., iOm) = / (X, Ki (X),..., (X),ccj+i,..., Wm) we get from ( |A.29| ) that 

dRj {pj+i {x, t) ) _ dRj {pj+i {x, t)) 


dx 


dp 


-Qj {Pj+i{xR),Uj+i{xR),.. .,u„^{xR)) 


< Rj{pj+i{x,t))+ aio+ 3 i^j_i) { \ui{x,t)\j , xe[Dj,Dj+i]. (A.30) 

\*=i+i 


Therefore, employing (A.28) we get that 


Rj{pj+i{x,t)) < e-^^+ia 9 + 3 Q_i) (llpj(f)lloo)+ ^i+ie^^+'aio+ 30 -i) ( ^ 


lUj 


. (A.31) 


\i=j+l 


Using (A.28), (A.27), and the properties of class JC functions the proof is completed. 


Proof of Lemma 

The proof of this lemma employs similar arguments to the proof of Lemma with the difference that 
one uses the ODEs in x for tt*, f = 1,..., m together with Assumption We again prove this lemma by 


induction. We first prove ( [63| ). Using the fact that tti satisfies the initial value problem ( |A.2[ ), ( |A.3| ), we 
get under Assumption and the definition of input-to-state stability (see, for example, [|5^ 1 that there 
exist a class ICC function and a class /Coo function 72 such that 

ki(x,f)| </33(|X(f)|,x)+72 f jap |M;i(x,f)| j V xg[0,L>i], (A.32) 


and hence, we arrive at (63) with pi(s) = jd^, (s, 0) + 72 (s). Assume next that for some j it holds that 


|7rj(f)||oo < Pi 


(A.33) 


Using the fact that satisfies the initial value problem (A.5), (A. 6 ) we get under Assumption that 

X G [Dj,Dj^i], (A.34) 

and hence, with (|A.33|) we arrive at ||7rj+i(f)||oo < Pj+i with Pj+i{s) = fds {pj, 0) + 72 (s). 


\7ij+i{x,t)\ < fd3{\7ijiDj)\,x - Dj)+-f 2 \ sup V |u;i(x,f)| 

\Dj<y<x 


Proof of Lemma 


We prove first ( 66 ). Using the fact that Ki, i = 1,... ,m, are locally Lipschitz with Ki(0) = 0, / 
1 ,..., m, there exist class /Coo functions a*, z = 1 ,..., m, such that 


|Ki(X)| < a* (|X|), i = l,...,m, 


(A.35) 
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for all X G Hence, using 


and relations (60)-(62) from Lemma we get estimate 
with p(s) = s + (j2j=T Similarly, using (52 >-(54) we get estimate (67) by employing 

Lemma 1^ with p(s) = s + (EjZT 
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